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Abstract. - We study spin precession due to Rashba spin splitting of electrons and holes in 
semiconductor quantum wells. Based on a simple analytical expression that we derive for the 
current modulation in a broad class of experimental situations of ferromagnet/nonmagnetic 
semiconductor/ferromagnet hybrid structures, we conclude that the Datta-Das spin transistor 
(i) is feasible with holes and (ii) its functionality is not affected by integration over injection 
angles. The current modulation shows a universal oscillation period, irrespective of the different 
forms of the Rashba Hamiltonian for electrons and holes. The analytic formulas approximate 
extremely well exact numerical calculations of a more elaborate Kohn-Luttinger model. 



Transport effects based on coherent manipulation of the spin degree of freedom in low- 
dimensional semiconductors are currently attracting a lot of attention [1]. These studies, 
enabled by recent progress in nanofabrication technology to create high-quality samples, are 
motivated by both their interesting fundamental physics and their potential for future device 
applications [2] . A lot of progress in the field has been stimulated by the exploitation of spin 
precession due to Rashba spin-orbit (SO) coupling in 2D systems both for electrons [3-5] and 
for holes [6]. A prominent example is the spin-controlled field-effect transistor (spin FET) 
introduced by Datta and Das [7] , followed by more recent proposals for novel devices utilizing 
Rashba SO coupling [8]. Both in the original [7] and most subsequent [9-14] works, a quasi- 
onedimensional (ID) confinement was considered essential for proper spin-FET action. Spin 
precession in truly 2D electron systems was studied numerically in a number of works [15,16]. 
On the other hand, ever-present spin relaxation will reduce the spin polarization of currents, 
preventing the realization of gate-controlled modulation. Such processes arise, e.g., from 
magnetic impurities but most importantly from elastic impurity scattering that randomizes 
the direction of the effective Rashba field. Stronger spin-orbit coupling and band mixing 
phenomena imply a shorter spin relaxation time for the holes respect to the electrons, that 
can be compensated by shorter precession length. A nice proposal, which exploits tunability 
of Rashba SO coupling, to overcome the detrimental effects due to scattering processes is 
presented in the last paper of Ref. [8] . 
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Fig. 1 - Current modulation in a p-type 2D ferromagnet/nonmagnetic-semiconductor/ferromagnet 
hybrid system computed numerically (solid line) for a realistic system and compared with our analytic 
expression (dashed line). Here L is the width of the nonmagnetic region, and Lso the spin-precession 
length for holes due to the Rashba effect. The magnetizations of the two ferromagnetic regions are 
parallel and point in a direction perpendicular to the 2D plane of the hybrid system. For the numerical 
simulation we have used parameter values for a realistic GaMnAs/GaAs/GaMnAs 2D heterostructure. 
The hole densities are of order 10 16 m" 2 , L$ « 300 nm, and the exchange field is 150 meV. There 
is a Fermi wave vector mismatch fey /&f,o = 1-33; and the Fermi wave vector splitting due to SO 
coupling is Afc/fcp,o = 7.66 • 10 -3 , corresponding to f3 p {E z ) — 1.3 ■ 10 20 eV _1 sec -2 nm. 



The aim of this Letter is twofold. First, we provide a unified analytical description of spin 
precession for both electrons and holes in the case of ballistic transport regime, discussing com- 
mon universal features and retaining the 2D nature of the problem. This is the most realistic 
case for anticipated device applications. One of our main results is illustrated in Fig.n wncrc 
our approximate analytical expression for the current modulation in a 2D hole spin FET is 
compared with the full numerical result obtained by mode-matching within a more elaborate 
Kohn-Luttinger model. Except for very small values of the distance L between the ferromag- 
netic contacts, the agreement is excellent. We prove analytically that the oscillation period 
in the 2D setup is the same as in the quasi-lD structures considered before [7,9]. Second, 
we want to emphasize the utility of holes as carriers in a spin FET. The apparent disadvan- 
tage of a shorter spin life time as compared to electrons is off-set by shorter spin-precession 
lengths and, above all, by the possibility to use recently created [17] p-type ferromagnetic 
semiconductors as current injectors. Besides obvious advantages concerning integrability in 
current semiconductor technology, such an all-semiconductor spin FET would circumvent the 
problems [18] that prevent spin injection from metallic magnets into semiconductors in the 
absence of a large interface barrier. 

We now turn to explaining details of our calculations. A 2D hybrid system with two semi- 
infinite ferromagnetic contacts separated by a nonmagnetic 2D stripe of width L is considered. 
To be specific, the extension of both ferromagnetic and nonmagnetic parts is infinite in y 
direction, and the nonmagnetic stripe is bounded by interfaces to the ferromagnetic contacts 
at x = and x = L. Due to structural inversion asymmetry in the growth (z) direction, 
charge carriers are subject to a spin-orbit coupling of the Rashba type [3,5,6]. For electrons 
in the s-like conduction band, it reads (3 s (E z )(px a) ■ z, while for the p-like valence bands we 
have ftp (E z )(p x J) • £ [19, 20] . Here, (E z ) is related to the average of the electric field in growth 
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Fig. 2 - a) Plot of the functions F(-y) and G < -°' h - ) (7) computed numerically (solid lines) and by means 
of the approximated formulae given in the text (dashed line). The difference between analytical and 
numerical results is only barely visible for 7^0. b) Current density as a function of 7 = nL/L^ 1 ' 
for different magnetization directions for both the electron and hole cases. The current is computed 
using the approximated expressions for f (7) and G^ e ' h \-f). 



direction. In 2D hole quantum wells, the degeneracy of heavy-hole (HH) and light -hole (LH) 
valence bands at the zone center is lifted by the vertical confinement [21], and only the first 
HH subband (HHi) is populated for typical 2D sheet densities. In the following, we consider 
only the HHi subband unless specified otherwise. The effective Hamiltonians for Rashba SO 
coupling experienced by electrons and holes are [19,20] 

= f3 s (E z )ilp_a + -p + a-} (la) 
i?io h) = ME z )i [p 3 _a+-p 3 + a.] , (lb) 

where p± = p x ± ip y are linear combinations of momentum components, a± = [a x ± ia y )/2 
denote Pauli matrices, and /3 s ,h are material parameters. 

We assume identical ferromagnetic contacts having parallel magnetization M whose direc- 
tion (in 3D) is described by two polar angles (0, </>), where = corresponds to magnetization 
perpendicular to interface. As quasi-2D ferromagnets have either easy-axis or easy-plane 
anisotropy, the two relevant cases are that with magnetization perpendicular to the 2D plane 
(9 = 0) or in the plane (9 — 7r/2). Assuming that there is no reflection at the interfaces 
between the nonmagnetic semiconductor and the ferromagnetic contacts, we can calculate the 
transmission probability Tg^a) for an electron impinging on the first interface at an angle 
a. (Details and a discussion of the approximation are given further below.) The total current 
density is then J oc J^^ 2 ^8 ,4>{ a ) cos(a) da. We find its expression for arbitrary (9, (f) as 

J = J { cos 2 6F(i) + sin 2 9 sin 2 + F( 7 ) cos 2 <j) + G (e ' h) (7) cos(20) \, (2) 

where the functions F and G^ e ' h - ) are shown in the panel a) of Fig- El an d are defined explicitly 
below. Here 7 = 7rL/Lso' h \ where Lso = 2n/Ak is the spin precession length for electrons 
or holes, with Afc being the difference of Fermi wave vectors for the spin-split 2D subbands: 



^-StSSUw.) ' (3a) 



h \ 2 1 
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Here m( c,h ) is the effective mass for electrons/holes, and £f is the Fermi energy. Equation (|3b|) 
is valid when the Rashba spin splitting is small compared to the Fermi energy. Equation J5J is 
one of the central results of this Letter, and a few remarks on it are in order. Once we decide 
whether the carriers are electrons or holes and, hence, fix the symmetry of the band (s or p 
type), the current density has a universal behavior: it depends only on the ratio of the distance 
between the ferromagnetic contacts and the spin-precession length. Finite transparency of the 
interfaces and Fermi-velocity mismatch between the ferromagnetic and nonmagnetic materials 
will be shown below to lead only to small quantitative changes. The only difference in the 
functional form of the current for electrons and holes is in G' ' 11 ' (7) and thus appears only if 
the magnetization has an in-plane component. The current density Eq. shows oscillations 
in 7 with a period 7r, i.e., oscillations in L with a period Lit , that are a manifestation of spin 
precession in the nonmagnetic region. A few observations can be made regarding the behavior 
of the oscillatory part of the current density for different magnetization directions (see panel 
b) of Fig.[2J): i) The largest oscillations of the current (solid line in the plot) as a function of 
7 occur for magnetization perpendicular to the 2D system (9 = 0); ii) In this case electrons 
and holes behave exactly in the same way; iii) For in-plane magnetization perpendicular to 
the interface (9 = ir/2, <fi = 0), the electron system (dashed line in the plot) exhibits a slightly 
larger modulation of the conductance as compared to the hole system (dash-dotted line in 
the plot); iv) For in-plane magnetization parallel to the interface (9 = 7r/2, (j> = 7r /2) the 
conductance still shows oscillations as a function of 7 (this is in contrast to the case of normal 
incidence), which are larger for holes than electrons. The different functional form for the SO 
term for electron and holes [see Eq.lfl]l] does not affect the transmission probability for normal 
injection. Hence, the non universal features, when present, are due to the integration over the 
injection angle a. We can conclude that the oscillations of the conductance as a function of 
Lso' 1 ^ are not washed out by the 2D geometry of the system, i.e., by the necessity to integrate 
over the direction of the incoming electrons in the ferromagnet. This oscillatory behavior of 
the current is exploited by the spin FET. Hence, our results show analytically that in the case 
of ballistic transport regime the spin FET does not require a quasi-lD setup to work properly, 
and that the largest modulations of the conductance as a function of the spin precession length 
are obtained for magnetization perpendicular to the 2D system. Furthermore, Eq. @ proves 
the feasibility of a p-type spin-FET, with holes performing even better than electrons for 
certain magnetization directions. 

To illustrate the strategy employed to obtain the general formula Eq. |3J) without inflicting 
lengthy algebraic manipulations to the reader, it is sufficient to show the calculation for the 
case of a hole system and magnetization perpendicular to the plane (9 = 0). We start by 
calculating the probability for a plane wave in the left ferromagnet, impinging on the first 
interface at an angle a, to be transmitted in the second ferromagnet. We assume that only 
majority spins contribute to transport, hence the spin state of the incoming electron in the 
left ferromagnet is |+), where |+) is the spinor corresponding to spin in the +z direction 
(the magnetization direction) . We now need to write the wave function in the semiconductor. 
We notice first that the SO coupling Hamiltonian Eq. (|lbl) removes the degeneracy of HHi 
subbands. The eigenstates are still plane waves exp [ik(x cos a + ysina)] x±j with spinors 
X± = (l/v / 2~)[|+) =p i exp(i3a)|— )]. The spin-split dispersion relations associated to these 
states read e± = 2 ^<h) k 2 ± (3h(E ' z )h 3 'k 3 . As we study linear transport, we are interested in 
the states at the Fermi energy. In particular, the Fermi wave vectors for the e± bands are 
fcp,± = k F ^Ak/2 where fc Fi0 is the Fermi wave vector when no SO coupling is present. Due to 
translational invariance along the interface, the wave vector component parallel to the interface 
is conserved when going from the ferromagnet to the non magnetic semiconductor. A plane 
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wave in the ferromagnet with k = fci, (xcosa + ysina) gives rise to two transmitted waves 
with propagation directions defined by angles a±. This effect is similar to birefringence [10, 
22]: the bands e± have different Fermi wave vectors, resulting in two different propagating 
directions for the transmitted waves. In the limit of weak SO coupling, i.e., Ak/k^ <C 1, 

the two angles read a± = ao ± (l/2)(Afc//cF,o) tanao, and ao is defined by ky sin a = 
fcF,o sin ao- We are now in the position to write the transmitted state in the nonmagnetic strip: 
it simply reads c+ exp [ik-p, + (x cos a + + ysma + )] x+ + c_ exp [ikF,-(x cosa_ + ysina_)] x_. 
By assuming perfect transparency of the interface we can compute the coefficients c± simply 
by matching the wave functions in the ferromagnet and non magnetic semiconductor. At the 
other interface x = L, only the |+) component will be transmitted, hence the outgoing state 

in the right ferromagnet reads exp ikjP {x cos a + y sin a) cos[A/cL/(2 cos ao)] \+) . From that 
the transmission probability can be read off: 



cos 



7 



cos a 



(4) 



where we have used the relation Afc L/2 = 7, and the dependence on a is through ao via the 
relation fc F sin a = kpfiSinao. In a similar way we can obtain the transmission probabilities 
for the electron case and for arbitrary magnetization direction. We find that Eq. (0J is valid 
for both electrons and holes. The transmission probability for in-plane magnetization reads 



T 



tt/2,0 



(a) 



cos 



7 



cos ao 



+ sin [v e ^a 



7 



cos ao 



(5) 



where v c = 1 and i\ = 3. Finally, we can write the transmission for arbitrary magnetization 
direction as 

,2 , _^_2 . 



Tg^(a)=cos 6>T Oi + sin 0T w / 2 ,, 



(6) 



Equations |0][|)J) cease to be valid once one of the transmitted states in the nonmagnetic strip 
becomes evanescent and is totally reflected. This condition defines the critical angles a C; ±, 
that in the limit of weak SO coupling read a, 



o,± = Ttrr T » TCff = «c At this point a few 

rc-p fcp 

remarks on the approximate analytical treatment are in order: i) In this calculation scheme the 
spin precession in the nonmagnetic strip and the spin selecting properties of the ferromagnets 
(they act as polarizer and analyzer) are fully taken into account, while the interference effects 
arising from multiple reflection between the interfaces are neglected. This is in the same 
spirit of the calculation of Datta and Das [7] and of all the quantum-mechanical Gedanken 
experiments involving polarized photons and polarizers. Indeed, multiple reflection introduces 
a modulation of the transmission coefficients oc cos 2 {Lk-p / cos a) . ii) This fast oscillation is 
washed out by the integration over the injection angle. This explains the remarkable agreement 
with the full quantum-mechanical calculation, shown in Fig. ^ 

The current density perpendicular to the interface is proportional to J"^ Tg^(a) cos a da oc 

/-7r/2 Te,<f>{oto) cosao dao, where Tg^(ao) = Tg t< f,(a(ao)). Writing this integral with Tg^(a) 
given in Eqs. 10]©, we obtain Eq. (J2J, where ^(7) and G^ c ' h ^ are 



^(7) = 2 J-f cos a cos2 (sh) da > 
G(e,h)( 7 ) = 1 cos a sin 2 (^, h a) sin 2 (^) da. 



(7) 
(8) 
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We have obtained an approximate analytical solution of Eqs. (Q and iJHJ that is valid unless 
7 is very small: 



F( 7 ) 



^{ (47 - 1)S ( 2 ^) + (-4 7+1 , C ( 2 ^l + 



:2a /- [cos(2 7 ) + sin(2 7 )] + 1 + — -= 

TT 7TV2 



(9) 



G(°' h )( 7 ) = V2^|p (c ' h '( 7 )C^2^+p( c ' h )(- 7 )S^2^- 

Q (c ' h) ( 7 )cos(2 7 ) + Q( c ' h) (-7)sin(2 7 )l +Af( c ^( 7 ) L ( 1.0) 



where C and S" are cosine and sine Fresnel integrals [23], respectively, and 

n(h) , . 61 7 16 6 4 5 4 11 3 7 2 9 

PW( 7 ) = 7 7 7 b + - 7 5 - 7 4 H Y Y H 7 - 

w; 315 ' 45 ' 3 ' ' 12 ' 16 ' 64 ' 



13 

256 ' 



)(e) 
)(h) 



1,1 1 

-7 + -7 H 

3 ' 6 ' 16 



32 



16 



218 4 34 3 11 2 7 



13 



-7 H 7 H 7 H 7 H 7 H 7 + 

315 ' 105 ' 315 105 ' 24 ' 48 ' 128 



M^( 7 ) = - 



7 



32 6^x72 ' 



nrth), \ 16 6 4 7 2 !7V2 13 
M(h, W = 45 7 +7 + I6 7 + 140^ 256 • 



(11a) 
(lib) 
(11c) 
(lid) 
(lie) 

(llf) 



In Fig. [3 we compare the approximate results Eqs. JHJ and (|10|) with the integrals of Eqs. Q 
and (jHJ- The approximate formulae work extremely well: the largest errors are for 7 = 0, 
and hence not in the region of physical interest, and the errors decrease monotonically with 
increasing 7. The relative errors for 7 = ir, i.e., after a full precession, are of order 10~ 3 [24]. 

We have derived the above results for transmission probabilities [Eqs. ©] and the current 
density [Eq. J2J] making use of several approximations, such as perfectly transparent interfaces 
and, in the hole case, the truncation to the lowest heavy-hole subband. We now test the 
accuracy of our approximate results by comparing them to a more realistic model for the 
hole system which is the more complicated one. Both the magnetic and nonmagnetic parts 
of the hybrid system are modeled by a 4 x 4 Kohn-Luttinger Hamiltonian with HH and 
LH bands. The confinement in the growth direction is treated by means of the envelope- 
function approximation. The ferromagnetic contacts are described within the Stoner approach 
by an exchange field. Calculation of transmission and reflection coefficients by means of a 
full quantum mechanical mode matching at the two interfaces and integration over angle 
of incidence yields the current density shown in Fig. ^ as a solid line. Excellent agreement 
between our approximate analytical results and the exact numerical curve is apparent (note 
that for the latter a Fermi wave vector mismatch across the interface is present.) The only 
effect of nonidcal interfaces on the current modulation is a renormalization of its amplitude 
Jo- 

In conclusions, we have studied spin precession in electron and hole system, finding a 
universal expression for the current in a spin FET geometry. The full 2D nature of the 
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problem is retained. Our results show that the current through the system is determined only 
by the common magnetization direction in the ferromagnets and by the ratio of the distance 
between the ferromagnetic contacts and the spin-precession length. 
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